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Abstract 

We study the symmetric Anderson-Holstein (AH) model at zero temperature with Wilson's 
numerical renormalization group (NRG) technique to study the interplay between the electron- 
electron and electron-phonon interactions. An improved method for calculating the phonon prop- 
agator using the NRG technique is presented, which turns out to be more accurate and reliable 
than the previous works in that it calculates the phonon renormalization explicitly and satisfies 
the boson sum rule better. The method is applied to calculate the renormalized phonon propaga- 
tors along with the electron propagators as the onsite Coulomb repulsion U and electron-phonon 
coupling constant g are varied. As g is increased, the phonon mode is successively renormalized, 
and for g > g co crosses over to the regime where the mode splits into two components, one of which 
approaches back to the bare frequency and the other develops into a soft mode. The initial renor- 
malization of the phonon mode, as g is increased from 0, depends on U and the hybridization A; it 
gets softened (hardened) for U > (<)U S (A). Correlated with the emergence of the soft mode is the 
central peak of the electron spectral function severely suppressed. These NRG calculations will be 
compared with the standard Green's function results for the weak coupling regime to understand 
the phonon renormalization and soft mode. 
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I. INTRODUCTION 



An important problem in interacting many particle systems is understanding the inter- 
play between the electron-electron and electron-phonon interactions. Such interplay will be 
important to the problems such as the effects of phonons in the strongly correlated elec- 
tron systems and the effects of electron correlations in the electron-phonon coupled systems. 
They include the cuprates, organic superconductors, Kondo systems, fullerenes, and the 
transition metal oxides, among others. Although some perturbative approaches have been 
proposed to understand the interplay M , there has not been much progress in this field par- 
tially because there are no reliable calculation schemes applicable in strong coupling regime 
where standard perturbation methods break down. Some non-perturbative techniques have 
been developed for strongly correlated impurity model in the context of the Kondo prob- 
lem, but these techniques can not be applied to the more interesting correlated models on 
a lattice. With the recent advent of the dynamical mean-field theory (DMFT), however, 
those techniques can be and have been successfully applied to a wide class of lattice mod- 
els. Within the DMFT, a lattice model is mapped onto an effective impurity model || [|. 
Among the various methods applied to solve the effective impurity model, Wilson's numer- 
ical renormalization group (NRG) is particularly powerful in that it is non-perturbative in 
nature so that it can cover the whole parameter space, and it can be applied to both zero 
and finite temperatures and to both the dynamical and static properties ||, EL BL [/]]. The 
NRG+DMFT technique has been successfully applied to the half-filled Hubbard model to 
study the metal-insulator transition at both T = and finite temperatures || . 

The prototype model for describing the interplay between electrons and phonons is the 
Hubbard-Holstein model and the Anderson-Holstein (AH) model which is the impurity ver- 
sion of the Hubbard-Holstein model. Recently, the AH model was studied at T = with the 
NRG technique by Hewson and Meyer ||11||. The Holstein model was also studied with the 



DMFT combined with the NRG technique [12]. These works reported electron and phonon 
spectral functions among others. In the present work, we focus on the phonon propagator 
and report an improved method to calculate the phonon spectral function using the NRG 
technique for the AH model at zero temperature. This is the necessary first step for applying 
the DMFT to the Hubbard-Holstein model. In section II, we present the improved method 
for calculating the phonon spectral function, which turns out to be more accurate and reli- 
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able than the previous methods in that it calculates the phonon renormalization explicitly 
and satisfies the sum rule better fTI|| . We will also present the standard Green's function 
analysis for the AH model, which provides a useful reference to understand the NRG results 
for the weak-coupling regime. Then, in section III, the results of these NRG calculations 
will be presented for the phonon and electron spectral functions, and they will be compared 
with the standard Green's function results to understand the phonon renormalization and 
soft mode. Finally, section IV is for the summary and concluding remarks. 



II. FORMULATION 



The AH model is defined by 

U AH = E wla** + e f £ ftf a + Uf\hf\f i + lu (a^a + ±) (1) 

+ £ Vk {Hcka + clfa) +9(a^ + a) (E flU - l) , 

fc,cr V a J 

where a is the spin index, ef the impurity level energy, and Vf. is the hybridization matrix 
element of the impurity f a electron with the conduction electrons of the host metal described 
by the operators Ck a - The / electron has the onsite Coulomb repulsion U, and is linearly 
coupled with the Einstein phonon of frequency Uq with the coupling constant g at the 
impurity site. The coupling is taken to be proportional to (J2 a ft fa ~ 1); where the factor of 
— 1 renders the problem particle-hole symmetric for the symmetric AH model, that is, the 
electron spectral function which will be discussed below is an even function of to. The AH 
model of Eq. ([!]) can be rewritten as 

n AH = (c, + |) E flf* + | (E flU - l) + (2) 
+ E tkcl a c k<7 + E Vk (ftcka + c\j a ) + g (a) + a) I J^ftU ~ 1 J , 

k,a k,a \ cr / 

disregarding constant terms. We will focus on the symmetric case of e/ + U/2 = in the 
present work for simplicity. Asymmetric cases are also straightforward. 

The electron and phonon Green's functions on the impurity site are defined, respectively, 

by 

G^) = ((U,ft)U d(u) = ((a,a%, (3) 



where the correlator (( )) is defined as 

/CO 
dte™ «Oi,O a )) t , (4) 
-co 

((o 1 ,o 2 )) t = -iB{t)([o 1 (t),o 2 (o)] c )^-ie(t)([o 1 (o),o 2 (-t)] c ). 

Here, 6 is the step function, ( ) the thermodynamic average, and ( = + if both 0\ and 2 
are fermion operators and ( = — otherwise. The equation of motion (EOM) of the correlator 
is given by 

w«0i,02» w = ([Oi,0 2 ] f ) + <([Oi,7fl_,0 2 » w , (5) 



The /-electron Green's function can be written using the EOM [11, 13 



G-\uj) = uj - e, - A(w) - ^ (w) - E»( W ), (6) 
where the self-energy corrections are given by 

This form of defining G a (u) was proposed and its advantages over the direct NRG calculation 
°f ((/) P))u were discussed by Bulla et al. and Hewson and Meyer ]lT], [L3[ . 



Now, we consider the phonon propagator d{uj) defined by Eq. (|3]). There are several ways 
to calculate the phonon propagator: (a) We may calculate d{uS) = ((a, a))) u directly by NRG, 
or, (b) calculate 



where 



d Q 1 (u)) = w-u Q , n f = J2flfa -1, (9) 



or, (c) calculate 



(10) 

Eqs. (H) and (ITOf) may be derived using the EOM as follows. From the Eq. fl5|), we have 

(u-u Q )^a,a^)) u = 1 + fi'((n / ,a t )) w . (11) 



Eq. (|IID can be rearranged, by dividing both sides by ((a, cfi)) u , to Eq. @. Again from the 
EOM, we also obtain 

uj{{n f ,a^))u = uo((nf,a*)) u + </«n/,n/)) w . (12) 

Substituting this into Eq. flTT|) yields Eq. (0). 

The advantage of using Eq. (^) is that it is written in the form of Dyson equation and 
shows the renormalization explicitly. The form of Eq. (pi]) was employed previously [|ll], 12| . 



It will be extremely difficult, however, to obtain the renormalized spectrum near u ~ uj 
from Eq. (|To|) because of the strong singularity of the [rf (^)] 2 - We therefore use Eq. (H) 
in the present analysis. This argument is made more quantitative by the sum rule. The 
important point is that the d(u) must satisfy the boson sum rule: 

1 r°° 

du Ira d{u) = 1. (13) 



7T 



oo 



The three ways of calculating the phonon propagator above will be discussed in the following 
section. It will be shown that the method (b) satisfies the sum rule of Eq. (|i~3| ) much better 
than the other methods. 

Alternatively, one may use the phonon Green's function defined by 



D(u) = ((a + a!,a + af)) u . (14) 



We then have 



D-H-) = D^)-g , (15) 

{{a + a\ a + al)) u 

D(u) = A)H + \gDo(u)] 2 ((n f ,n f )) u , 

which correspond to Eqs. (§) and (^0|), respectively, where Dq(uj) = 2ujq/{uj 2 — Uq). There 
are essentially no differences between the phonon Green's functions d(u) and D(u). The 
advantage of d(uj) over D{uj) is that there exists the useful sum rule of Eq. ( fL3|) the d{uj) 
satisfies, while the sum rule the D(u) satisfies 

- duImD(uj)n(uj) = ((a f + a)(a f + a)) , (16) 



TC J-oo 

where n(u) is the Bose function, is less useful. 

Before we present the detailed NRG results for the AH model, we consider the standard 
Green's function analyses of the model for the phonon and electron propagators. It will 
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FIG. 1: The Feynman diagram for the Anderson-Holstein model. The single (double) solid and 
wiggle lines represent, respectively, the bare (renormalized) electron and phonon propagators. The 
single dashed line refers to the Hubbard interaction, and the double dashed line to the effective 
interaction between two electrons of spin a and (3. 

provide a useful reference to understand the NRG results for the weak-coupling regime. 
Within the perturbative Green's function scheme, the renormalized electron and phonon 
propagators for U = shown in Fig. 1(a) and (b), given by 

G- 1 M = G 1 M-E( W ), D-\u) = D^(uj)-g 2 x(uj), (17) 



where Go is the bare electron Green's function given by 

1 

uj + iA 



G (uj) = ^—, A = Y,*\V k \ 2 5(u J -e k ), (18) 



k 

and the self-energy £ and the polarization % are given by 

= -L [°° de r dQ f ^\ + n (~ Q \ lmG(e)ImV(n), (19) 
it J-oo J-oo e + il — uj — id 

X {u) = - 2 r de r de' f [ e) ~ f , ie ' ] - ImG{e)ImG{e'). 

7T Z J-oo J-oo uj + e — e + to 

A of Eq. (|18D is the minus imaginary part of A defined in Eq. (^), A = — iA. V is the 
effective interaction between two electrons which is reduced to V = g 2 D for U = as 
discussed below. We consider the case where A (a;) is a constant. At T = and up to the 
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2nd order in g, we obtain 

xH = 



1 



. ■ In I 1 — i— I . 



+ 



1 r 00 A 

7r Jo e 2 + A 2 Ve + a; + zA ' e — — iA 
2A 1 / 



(20) 



The real and imaginary parts are given, respectively, by 



Imx(w) = 
For uj <C A, it is reduced to 



A 



% 4A 2 + uj 2 



CO ' 



\ 4A 



UJ 



ln|1 + A^J-i7 tan "A 



2A 2 



1 



■nui 4A 2 + cu 2 
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In 1 + — + — tan" 1 — 
V A 2 / A A 



5 a; 2 



Xl^J 



1 - 77^7 +^ 



(21) 



(22) 



ttA V 6 A 2 ' "A y 1 ' 
which will be used in discussing the soft phonon mode in the following section. 
When U ^ 0, the renormalized phonon propagator is given by 

This may be obtained as follows. The effective interaction between the /-electrons V a /3, 
where a and (3 are the spin indices, due to both electron-phonon and electron-electron 
interactions is determined, as shown in Fig. 1(c), by 



where 



a', /3' 



C ee V° = g*D , V°=U + V°. 



(24) 



(25) 



The solution to Eq. p3) is given by 



Vn = Vv 



V° + U(2V° + U)xk 



[1 - (2V° + C/)(xtt + Xn)l [1 + ^ (Xtt " Xu)] ' 

v° + u - u(2V° + u) Xn 



n iT [i - (2V° + U){ xn + xu)] [i + ^(xtt - xu)] ' 

We put xtt = Xu = X/2 and xn = X|T = to obtain 



(26) 



yo + ^(yo + ^ 2 )x 



[i-(y° + c//2) x ] [i + c/x/2]' 

[1 _ (yo + [i + c/ x /2] 
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(27) 



The renormalized electron-phonon interaction is then given by |10 



U 2 x/2 U 

Vel-ph = V n - - — — = V U - 



l-U 2 X 2 /4: n l-f/ 2 x 2 /4 
V° 

(28) 



- [l-(2yo + C /) x /2][i_[/ x /2]- 
We put V° = g 2 D Q to obtain 

= h jf ( W91 2 ^M. (29) 
[1 - ?7x(w)/2] 

where D{u) and x(u;) are given, respectively, by Eqs. ( [2"3"D and fl2Tj|). We will use the 
renormalized phonon propagator of Eq. ( |23| ) in the following section to understand the 
phonon renormalization and emergence of the soft phonon mode. 

III. SPECTRAL FUNCTIONS AND SUM RULE 

The renormalized phonon Green's function d(u) was obtained by the method (b) of 
Eq. (|J) of the previous section by calculating ((nf,^))^ and ((a, a'))^ directly with the 



NRG technique. Some technical details of the NRG scheme can be found in Ref. [0, [13 
Briefly, the conduction electron band of the AH model, which was taken from —1 to 1, is 
logarithmically discretized and mapped onto a semi-infinite fictitious chain coupled to the 
impurity site at one end such that the hopping amplitude along the chain falls off expontially 
(~ A~ n / 2 ). The many body eigen-states and energies are calculated iteratively starting from 
the impurity site and adding an extra site along the chain at each step. After a few steps, 
the number of states has to be truncated to the lowest N st ones. In the present work, we 
used the discretization parameter A = 2.0, kept N st = 600 states in each iteration, and 100 
fictitious sites along the chain. The lowest 40 phonon states are kept for the AH model with 
the bare phonon frequency ujq = 0.05. The energy unit is the half the conduction bandwidth. 

We now compare the three ways of calculating the d(u) by the methods (a), (b), and 
(c), given, respectively, by Eqs. (|3]), (§), and (|10D explained in the previous section. In 
Fig. 2, we show —^Imd{uj) vs. u calculated by the above three methods for U = 0.1 
and A = uj /tt w 0.0159, for the electron-phonon coupling constants g = 0.01, 0.03, 0.05, 
0.07, and 0.09. The NRG calculations yield the discrete frequencies; the resulting spectral 
functions are discrete with the 5-peaks getting closer as uj — > 0. The 5-peaks are broadened 
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FIG. 2: The phonon spectral function —^Im d(u>) vs. u> calculated by the three methods for U = 0.1 
and A = uoq/tt ~ 0.0159. Figure (a), (b), and (c) are obtained by the method (a) the direct NRG 
calculation of ((a, cfi))^, (b) the present calculation of Eq. @, and (c) Hewson and Meyer of Eq. 
(|l0|) , respectively. 



via 



8(u - u n ) 



-b 2 /A 



exp 



(In a; — liau r , 
b 2 



(30) 



This broadening function is a Gaussian on a logarithmic scale with the width b. b = 0.4 was 
used in the present study. Figures (a), (b), and (c) are the results obtained by the methods 
(a), (b), and (c), respectively. For the low frequency, they show the same behavior, but for 
uo ~ uq, they show deviations from each other. In the method (a) the feature around the 
bare frequency ujq is almost washed out, while in the method (c) the peak is always tied to 
the bare frequency because of the singularity of the [g? (u;)] 2 . In the method (b), the phonon 
frequency is successively renormalized as one may expect from the perturbation theory. As 
the electron-phonon coupling constant g is increased, the phonon mode softens until the 
soft mode begins to emerge around g ~ g co , which will be discussed below. As g is further 
increased above g co , the mode splits into two components, one of which hardens back to the 
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g 


Direct NRG 


Improved method 


0.01 


1.000046 


0.999922 


0.03 


0.999212 


0.999890 


0.05 


0.993977 


0.999863 


0.07 


0.911200 


0.991001 


0.09 


0.992165 


0.999840 



TABLE I: The boson sum rule of Eq. (|i"3|). The second and third columns show, respectively, the 
integral — — du)Imd{u)) computed by the direct NRG calculation of (( (2, d ))t^ cincl th.6 improved, 
method of Eq. (§). 

bare frequency and the other develops into the soft mode as can be seen for g = 0.07 and 
0.09 of the figure (b). For g = 0.09, the soft mode can hardly be seen because the width of 
the soft mode is very narrow in this energy scale. 

In the table 1, we compare the integrals — - duilmd{uj) from the method (a) and (b) 
to check the sum rule of Eq. ([13|) . The method (b) satisfies the sum rule better than (a), 
and the method (c) doesn't even converge because of the strong singularity of the [g? (u;)] 2 
in Eq. (|T0|). Because of the singularity, it will be extremely difficult to obtain the phonon 
renormalization around ui ~ u> from the method (c). The Eq. (0) of method (b) is written 
in the form of Dyson equation where (l/g)((rif, a^^/^a, cr}^ is the proper polarization, 
which gives the phonon renormalization as g is increased. As g is further increased, the 
soft mode of u = emerges, which may be understood in terms of the Green's function 
perturbation scheme discussed below. The present method (b) of Eq. (H) calculated the 
phonon renormalization explicitly and satisfies the sum rule better than the other methods 
and, therefore, give more reliable and satisfactory description of the phonon renormalization 
of the AH model. 

The renormalized phonon frequency can be found from D~ 1 (oj) = 0. To investigate the 
emergence of the soft phonon mode, we substitute the expression Eq. (f22|) of x( u ) f° r ^ <^ A 
into Eq. (^) to obtain 

-^°-2jy=°- < 3i) 



47r^ 2 cJo(10vrA - U) 
3A(2ttA + Uf 
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FIG. 3: The electron spectral function —^ImG(uj) vs. ui for U = 0.1 and A = ujq/tt corresponding 
to the figure 2. 



The cross-over value g co is then given by 



9co = -^uj (27tA + U) 



(32) 



which gives g co = 0.05 for U — 0.1 and A = lu /it. We therefore expect the soft mode 
emerges for g > g co which agrees fairly well with the NRG calculations of Fig. 2(b). 

In Fig. 3, we show the electron spectral function A(uo) = —^ImG(uo) vs. uo for the same 
parameters as the Fig. 2. For small g such that U > 7rA, where U = U — 2E P is the effective 
onsite repulsion and E p = g 2 /ujQ is the polaron binding energy, the system is in the Kondo 
regime. The peak at u = is then the Kondo resonance peak whose width is given by 



and the height is 



27V~(C?A^ exp(--^-), 



A(u = 0) = -L, 



(33) 



(34) 



in agreement with the figure 3. The shoulders at uj ±U /2 xs ±0.05 correspond to the 
charge fluctuations at the / electron impurity site. As g is increased, U is decreased, and the 
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width of the Kondo peak increases as can be seen from the figure. The effective U becomes 



negative for g > g_ = yUio /2. g~ — 0.05 for U — 0.1 and lo = 0.05. The electron spectral 
function in this case may be understood in terms of the small polaron picture as discussed 
by Hewson and Meyer JD]. The electron Green's function for small hybridization limit of 
A — > is given by 



oo \n r 



n=0 n - 



((1 -n f>a ){l -n f _ a )) ((1 - 



oo — €f — nujQ oo — €f — U — uujq 

(35) 



(n/A 1 ~ n f,-o)) + ( n f,<T n f-*) 



oo — €f + rihjQ oo — €f — U + rwoQ. 

where A = g 1 '/oofi, (w/ )(7 ) = (fafa), an d tf = tf + E p . Note that 6f = 6f + E p instead 
of €f = €f — E p because the coupling between the / electrons and the phonons is taken 
to be proportional to (j2 a f\j a — l) instead of \ f^faj ■ The set of delta functions at 
oo = €f + nuo, €/ — nuo, e/ + U — nuo and e/ + U — nuo for integer n, determined by the 
denominators of Eq. (|35"|), will be broadened for non-zero A. The most probable phonon 
occupation number is (n) = E p /u = A. For integer n's, the peaks are expected at 

CO — €f ± TlLOo — €f + E p , €f + E p ± LOq, ■ ■ ■ , 

a; = 6f + U ± nt^o = tf + U — E p , ej + U — E p ± u , • • • , (36) 

which roughly tracks the four high energy peaks of the spectral function A(oo) as g 
is increased. For g = 0.07 for instance, E p = 0.098, and we expect the peaks at 
uo = ±0.048, ± 0.098, and so forth, which roughly agree with the detailed NRG results 
shown in the figure 3. The central peak for the negative U may be understood in terms 
of the Kondo-like resonance of the negative U model, which arises because of the charge 
fluctuations, instead of the spin fluctuations in the Kondo regime, between the zero and 
double occupation states at the impurity site. When g is further increased above g co , the 
soft phonon mode emerges. The effective \U\, which equals to 2g 2 /uoq — U for the Einstein 
phonon may be generalized to 

\U\ = 2# 2 (^) ~U, where = - J°° dco-Imd(oo)-. (37) 

When the soft mode emerges, the U can be very large and the width of the electron central 
peak given by Eq. ([33|) becomes very narrow. The central peak can hardly be seen for 
g = 0.07 and 0.09 of the Fig. 3 in this frequency scale. 
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FIG. 4: The phonon and electron spectral functions vs. uj for U = 0.2 and A = 0.016. 
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FIG. 5: The phonon and electron spectral functions vs. u for U = and A = 0.01. 
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We now present the results calculated with the method (b) for other parameter values. 
Fig. 4 is for U = 0.2 and A = 0.016 w uj /tt, and Fig. 5 is for U = and A = 0.01. They 
exhibit the similar behaviors as for the U = 0.1 and A = uj /tt of figures 2 and 3. The 
emergence of the soft phonon modes and suppression of the central peak are correlated for 
both cases and agree with the estimate of the perturbation theory of Eq. Q3"2]) as discussed 
above. g co = 0.0613 for U = 0.2 and A = 0.016 and g co = 0.0280 for U = and A = 0.01 
which are consistent with the NRG results of figures 4 and 5, respectively. The four high 
energy peaks for negative U are also consistent with the small polaron estimate of Eq. (j36l). 
For U = 0.2, shown in Fig. 4, e/ = —0.1, and E p = 0.128 for g = 0.08. We therefore expect 
to see the peaks in the electron spectral function at u ±(E p + e/), ± (E p + e/ + <^o) ; • • • ~ 
±0.028, ±0.078, and so forth, which agree with Fig. 4. For U — 0, shown in Fig. 5, = 0, 
and E p = 0.032 and 0.05 for g = 0.04 and 0.05, respectively. The peaks are expected at 
u « ±0.032 and ±0.082 for g = 0.04, and at u « ±0.05 and ±0.1 for # = 0.05, which are 
again consistent with the electron spectral function A(u) shown in Fig. 5. The heights of 
the central peak of the electron spectral function A(u) are given by 1/(ttA) of Eq. (|3"1D for 
both cases. 

For U = 0, the phonon renormalization seems somewhat different from the U ^ cases 
of Figs. 2 and 4 in that the main peak around uq is slightly hardened as g is increased. This 
apparent difference can be straightforwardly understood from the Green's function pertur- 
bation results presented in the previous section. If one solves for the renormalized phonon 
frequency from D^iuo) = using Eqs. ( p3|) and (]20f) , one can see that the renormalized 
frequency gets softened or hardened depending on A and U . The value of U s above (below) 
which the renormalized phonon frequency, as g in increased from 0, is softened (hardened) 
is simply given by Re[x(uJo)/(l — (U/2)x(wo))] = from Eq. (0), which yields 

U (A) - 2 RexM - 2 -^l ^ ^ + j) " 4 (j) tan " f38) 
4 J IXNP" A ln 2 (l + ^)+4tan-(^) ' 1 j 

For tuo = 0.05, one finds that £7, = 0.0636 for A = uj /tt, and 0.1868 for A = 0.01. 
We therefore expect that the phonon renormalization will get softened for U = 0.1 and 
A = ujo/iy, and U = 0.2 and A = 0.016 ~ c^o/vr, corresponding, respectively, to Fig. 2(b) 
and Fig. 4, but hardened for U = and A = 0.01 corresponding to Fig. 5, in agreement 
with the NRG results presented in those figures. 

Note, however, that the phonon spectral function of U = has a long tail towards ui = 
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compared with the U cases and, therefore, the average phonon frequency gets softened 
as g is increased. When g > g co , the soft mode emerges in addition to the peak around 
u , and the average phonon frequency hardens back to the bare frequency of u> as g is 
further increased above g co like the U ^ cases. Also noticeable is that the hardening or 
softening of the renormalized phonon frequency for U = does depend on A. From Eq. 

we expect that the renormalized frequency will get hardened (softened) for A < (>) 
0.410 u = 0.0205. The NRG results presented in Fig. 5 for U = and A = 0.01 are 
consistent with the estimate. 



IV. SUMMARY AND CONCLUDING REMARKS 

In the present paper, we have presented Wilson's numerical renormalization group cal- 
culations of the symmetric Anderson-Holstein model. The phonon and electron spectral 
functions were computed as the electron-electron interaction U and electron-phonon cou- 
pling g were varied. The electron spectral functions are in good agreement with the Hewson 
and Meyer results. On the other hand, an improved method was presented for calculating 
the phonon spectral functions and the results were compared with Hewson and Meyer's. 
Unlike the method (c) of Eq. fliCf ) where the phonon frequency is always tied to the bare 
frequency, the present calculations show that as g is increased, the phonon mode is renor- 
malized until the soft mode emerges around g w g co , and as g is further increased, the 
phonon mode splits into two components, one of which approaches back to the bare fre- 
quency and the other develops into the soft mode. We argue that the above behavior of the 
phonon renormalization is physically more sensible and is a correct picture in that (a) in 
the present method the phonon renormalization is explicitly computed, and (b) the present 
method satisfies the boson sum rule better than the previous methods as listed in Table 1. 
These NRG calculations were also compared with the standard Green's function results for 
the weak coupling regime to understand the phonon renormalization and soft mode. 

Now that the calculation of the dynamical properties of the impurity Anderson-Holstein 
model is successfully carried out, we plan to extend these calculations to the lattice model 
of the Hubbard-Holstein model within the framework of the dynamical mean-field theory. 
These calculations will be relevant to high-T c superconductors where recent photoemission 



experiments give evidence of strong electron-phonon coupling JTj] and to the manganites [15 
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Extension to the finite temperatures for the impurity model and for the lattice model should 
also be fruitful. 

We would like to thank Hyun Cheol Lee, Tae Suk Kim, Ralf Bulla, Dietrich Meyer, and 
Alex Hewson for helpful comments and discussions. This work was supported by the Korea 
Science & Engineering Foundation (KOSEF) through grant No. R01-1999-000-00031-0 and 
Center for Nanotubes and Nanoscale Composites (CNNC), and by the Ministry of Education 
through BK21 SNU-SKKU Physics program. 



[1] A. Deppeler and A. J. Millis, Phys. Rev. 65, 100301 (2002). 

[2] W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 (1989); ibid, 1066 (1989). 

[3] A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 

[4] H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. Rev. B 21, 1003 (1980). 

[5] H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. Rev. B 21, 1044 (1980). 

[6] K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

[7] R. Bulla, Adv. Solid State Phys. 40, 169 (2000). 

[8] R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

[9] R. Bulla, T. A. Costi, and D. Vollhardt, Phys. Rev. B 64, 045103 (2001). 

[10] G. D. Mahan, Many-Particle Physics, 3rd ed. (Plenum, New York, 2000). 

[11] A. C. Hewson and D. Meyer, J. Phys.: Condens. Matter 14, 427 (2002). 

[12] D. Meyer, A. C. Hewson, and R. Bulla, Phys. Rev. Lett. 89, 196401 (2002). 

[13] R. Bulla, A. C. Hewson, and Th. Pruschke, J. Phys.: Condens. Matter 10, 8365 (1998). 

[14] A. Lanzara, P. V. Bogdanov, X. J. Zhou, S. A. Keller, D. L. Feng, E. D. Lu, T. Yoshida, 

H. Eisaki, A. Fujimori, K. Kishio, J. -I. Shimoyama, T. Nodaand S. Uchida, Z. Hussain, and 

Z.-X. Shen, Nature 412, 510 (2001). 

[15] A. P. Ramirez, J. Phys.: Condens. Matter 9, 8171 (1997). 

[16] A. J. Millis, P. B. Littlewood, and B. I. Shraiman, Phys. Rev. Lett. 74, 5144 (1995). 



18 



